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We describe a method for initializing an ensemble of qubits in a pure ground state by applying 
collective cavity cooling techniques in the presence of local dephasing noise on each qubit. To solve 
the dynamics of the ensemble system we introduce a method for dissipative perturbation theory that 
applies average Hamiltonian theory in an imaginary-time dissipative interaction frame to find an 
average effective dissipator for the system dynamics. We use SU(4) algebra generators to analytically 
solve the first order perturbation for an arbitrary number of qubits in the ensemble. We find that to 
first order the effective dissipator describes local Ti thermal relaxation to the ground state of each 
qubit in the ensemble at a rate equal to the collective cavity cooling dissipation rate. The proposed 
technique should permit the parallel initialization of high purity states in large ensemble quantum 
systems based on solid-state spins. 


I. INTRODUCTION 

A fundamental challenge to implementing quantum in¬ 
formation processing on a physical device is the ability to 
rapidly and repeatable initializing the quantum system in 
a high purity state. Depending on the physical architec¬ 
ture of a device, quantum state initialization is typically 
done using methods such as strong projective measure¬ 
ments or filters [1], optical pumping and atomic transi¬ 
tions [2, 3], laser and microwave cooling [4-7], dynamical 
nuclear polarization [8, 9], algorithmic cooling [10, 11], 
and dissipative state engineering [12, 13]. As the num¬ 
ber of qubits in quantum devices increases initialization 
methods that can be implemented in parallel across many 
qubits are necessary to enable scalability. In the case 
of spin ensembles thermal relaxation processes naturally 
initialize all spins in parallel and can be used for state 
preparation, however to achieve high purity state at ther¬ 
mal equilibrium very low temperatures and strong mag¬ 
netic fields are needed. In the case of solid-state electron 
spin systems this can be a limiting factor as the thermal 
relaxation time (Tj) becomes very long at low tempera¬ 
tures [8]. 

In recent work by the authors and others it was sug¬ 
gested that cavity cooling techniques could enable the 
collective removal of entropy from an ensemble spin sys¬ 
tem [14, 15], and this has recently been experimentally 
demonstrated [16]. Our proposal in [15] relied on en¬ 
gineering a Tavis-Cummings (TC) exchange interaction 
between the spin ensemble and the side band of a high-Q 
resonator that was actively cooled to a low temperature 
thermal state. This is of particular interest as there have 
been several recent designs and proposals for the read¬ 
out and control of electron spin ensembles using high-Q 
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superconducting resonators [17-20]. Due to the identical 
particles in the ensemble the dynamics preserve a global 
SU(2) symmetry and the state space has a block diago¬ 
nal structure with each subspace block corresponding to 
an irreducible representation of SU(2) that behaves as an 
effective spin-J particle [21, 22]. The effective interaction 
experienced by the spin ensemble was dissipative cooling 
of each subspace to its respective low temperatures state. 
The largest effective spin subspace has J = N /2 and is 
called the totally symmetric subspace, or Dicke subspace, 
and contains all states which are permutation invariant 
across all N spins [23]. The cavity cooling dissipator de¬ 
rived in [15] preserves the global SU(2) symmetry, and 
so does not couple the subspaces. In this case, if cooling 
from a maximally mixed initial state the initial popula¬ 
tion in each subspace will be trapped in that subspace’s 
ground state, and the final state will be a mixed state sum 
of these ground states. The open problem is how to solve 
the dynamics of cavity cooling an ensemble quantum sys¬ 
tem while including an interaction which couples between 
these subspaces to enable cooling to the true ground state 
of the ensemble. The challenge is that once such an in¬ 
teraction is included, in principle the full Hilbert space 
for Wspins must be included which has dimension 4^. 

In the present paper we analyze a method for cavity 
cooling to the ground state of the ensemble by introduc¬ 
ing a local dephasing ( T 2 ) noise process on each spin in 
the ensemble. This acts to distinguish each spin and thus 
breaks the collective symmetry of the ensemble. In prac¬ 
tice dephasing is a phenomenological description of noise 
in solid-state spin systems can be caused by inhomoge¬ 
neous broadening, or spurious coupling to neighbouring 
spins or a local spin bath [8]. The effects of local dissi¬ 
pation on collective dynamics has been considered previ¬ 
ously, where it was shown to rapidly decohere coherent 
states in the Dicke subspace [22]. Decoherence of the 
Dicke subspace has also been studied for an inhomoge- 
neouly broadened ensemble of qubits coupled to a cav- 
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ity [24]. Unlike these previous studies we are able to use 
local dissipation as a resource for dissipative quantum 
state engineering to the ground state of the ensemble. 


Similar theoretical results to have been presented for 
cooling an ensemble of nuclear spins by coupling to the 
motion of a nanoscale mechanical resonator [14]. It was 
shown that the addition of a chemical shift to each spin 
was in principle sufficient to break the symmetry and 
achieve exponential relaxation to the ground state and 
this was numerically demonstrated for five spins. A sim¬ 
ilar result was found numerically in [25] where the relax¬ 
ation rates of the Dicke and non-Dicke subspaces were 
simulated for 10 inhomogeneously broadened qubits. Our 
work differs from the approaches in [14, 25] by includ¬ 
ing the symmetry breaking mechanism as a dissipative 
term leading to a Lindblad master equation that we can 
solve perturbatively to derive an analytic expression for 
the cooling dynamics in the regime where the first order 
perturbation term dominates the dynamics. In order to 
solve the master equation we develop a perturbation the¬ 
ory technique for dissipative evolution. This involves ap¬ 
plying the Magnus expansion [26, 27], or Average Hamil¬ 
tonian Theory [28-30], in an imaginary-time dissipative 
interaction frame to the superoperators describing evolu¬ 
tion. There is a long history of applying average Hamilto¬ 
nian theory to superoperators, called Average Liouvillian 
Theory [31-33], and also the related cumulant expansion 
approach for stochastic noise processes [34, 35]. In both 
of these cases the relevant interaction frame is defined by 
a Hamiltonian, which typically corresponds to a sequence 
of control pulses, and the average affect of a dissipative 
term in this frame is assessed. Our approach extends 
these formalisms by providing a procedure to apply these 
techniques in a purely dissipative (non-periodic) interac¬ 
tion frame. The utility of this method is that in a sys¬ 
tem with multiple decoherence mechanisms one can find 
the average effective dissipation of one mechanism in the 
presence of another. As we demonstrate, this may then 
be used for finding the equilibrium state of the system 
for dissipative state engineering applications. 


To analytically solve the lowest order term in the Mag¬ 
nus expansion we use recently introduced techniques for 
describing the local superoperators on spin ensembles us¬ 
ing SU(4) algebra generators [36, 37]. In this represen¬ 
tation the dynamics preserve a global SU(4) symmetry 
and any thermal state of the Hamiltonian with SU(2) 
symmetry will be entirely contained in the totally sym¬ 
metric subspace of SU(4). By considering the explicit 
matrix representation of the N qubit totally symmetric 
subspace of SU(4), this approach allows us to numerically 
simulate the reduced dynamics of the spin ensemble for 
up to N = 100 spins on a desktop computer. 


II. CAVITY COOLING WITH LOCAL 
DEPHASING. 

Consider an ensemble of N identical spin 1/2 particles 
interacting with a single mode cavity. Let a, be the 
cavity lowering and raising operators respectively, and 

N ( j) N 

J “ = £^r’ J ± = J2 a ± ) > (!) 

i=i 1=1 

be the collective spin operators for the spin ensemble, 
where ,a = x,y,zis the Pauli matrix for the j th spin- 
| system. If the spins are on-resonance with the cavity, 
or are driven to be on-resonance with a side-band of the 
cavity as was considered in [15], then the spin-cavity in¬ 
teraction is well described by the TC Hamiltonian 

#tc = g{J+a + T_a t ), (2) 

which is the N spin generalization of the Jaynes- 
Cummings interaction [38, 39]. The cavity will experi¬ 
ence photon-loss at a rate inversely proportional to Q, 
which may be described by the Lindblad dissipator 

V c = k(1 + n)D[a\ + KnD[af], (3) 

where 

D[a\p = apa f - ^{a'a,p}. (4) 

The effect of T> c is to reset the cavity to a thermal state 
p e q satisfying n = Tr[a'ap e q], at a rate k. In this case 
the evolution of the joint spin-cavity system is described 
by the Lindblad master equation 

^p(t) =-i[H TC ,p(t)\+V c p(t). (5) 

It was shown in [15] that by adiabatically eliminating 

the cavity in the Markovian regime (k gy/N), this 
interaction leads to an effective spin-dissipation of the 
form 

j t Ps{t) ='D cc p s {t), (6) 

where p s is the density matrix of the spin ensemble alone 
after adiabatic elimination of the cavity, and 

V cc = T(l + n)D[J-]+TnD[J+], (7) 

is the cavity cooling dissipator with rate T = 4 g 2 /n. For 
n = 0, the dissipator V cc cools each of SU(2) irrep sub¬ 
spaces of the collective spin ensemble to their respective 
ground states. However, since V cc preserves the SU(2) 
symmetry it does not couple the subspaces. Thus if the 
system is initially in a highly thermal state, cavity cool¬ 
ing dissipator alone is not sufficient to cool to the ground 
state of the system. 
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To enable cooling to the ground state we need to in¬ 
troduce an interaction that breaks the SU(2) symmetry 
of the ensemble, yet does not inhibit the cavity cooling 
dissipator. This may be achieved by a local dephasing 
term which acts identically on each spin in the ensemble. 
Dephasing, or X 2 dissipation, of a single spin-half sys¬ 
tem causes the off-diagonal density matrix elements of 
the spin to decay exponentially to zero. It is generated 
by the Lindblad dissipator 

D[a z /2]p = ^(a z pa f z - p). ( 8 ) 

The dissipator for identical local dephasing dissipators 
on each spin in the ensemble is then given by 

N 

V T 2 =^ 7 D[a«/2]. (9) 

3 =1 

This process results in a mixing of states across different 
spin-J subspaces that have the same J z value. Since cav¬ 
ity cooling drives each subspace to the lowest J z value 
state, a X 2 process will leak population trapped in the 
ground state of a spin-J subspace into the k th excited 
state of a spin-(J+£:) subspace, which will then be cooled 
to that subspaces ground state. In the ideal case this 
eventually leads to the ground state of the Dicke sub¬ 
space. To show this we must solving the the dynamics of 
the spin master equation 

j t Ps{t) = (V cc + V T 2 )p s {t). (10) 

A. Perturbative Solution 

We now introduce a new dissipative perturbation the¬ 
ory approach to solve Eq. (10). To begin we move into 
a dissipative interaction frame defined by the dephasing 
dissipator T>t 2 ■ In this interaction frame we have 

p s {t) =V cc (t)p s {t ), (11) 

where 

Ps(t) =e w ^p s (t), ( 12 ) 

V cc {t)=e tv ^V cc e~ tD ^. (13) 

Since this is a dissipative interaction frame, the time de¬ 
pendent terms in 'p s (t),'D cc (t) will be of the form e ±wt 
for real parameters w. Hence the terms e +ult diverge as t 
increases. By transforming into imaginary time (t 1 —> ir) 
we may convert this dissipative interaction frame into a 
periodic one leading to the periodic differential equation 

^(5(r) = «G(r)Q(r), (14) 

where 

(15) 


The explicit time dependence of the operator G(r) is 
given by 

G(r) =r(l+n)G_(r) + rnG + (r) (16) 

G± (t) = D[J ± ] + (e ±I7T - 1 )A± + (e Ti7T - 1 )B ± (17) 


where n is the average thermal photon number of the cav¬ 
ity, and A±, B± may be expressed in terms of SU(4) gen¬ 
erators (See Appendix A). Hence G(r) is periodic with 
period T = 27r/7, where 7 is the single spin dephasing 
rate. Eq. (14) has the formal solution 

Q(t) = Texp J dsG(s)j Q( 0 ), (18) 

where T is the time-ordering operator. If the dephasing 
rate is strong, then we can consider a perturbative expan¬ 
sion of Eq. (18) in terms of the Magnus expansion [26] 
to find an average time-independent dissipator over the 
period T: 



T exp (i dsG(s) j = exp ( iT D^ 


(19) 


fc=i 


Since the average dissipator is time-independent, we may 
transform back into real time to obtain an average de¬ 
scription of the system dynamics over the period T, which 
may then be used to compute the stroboscopic evolution 
over integer multiples of the period T: 


p s {nT) = exp 



Ps( 0). 


( 20 ) 


If the dephasing rate is greater than the collective cav¬ 
ity cooling rate we may make a secular approximation 
and only consider the lowest order term in the average 
dissipator Magnus expansion: 

D 1 = T- 1 [ dsG(s). (21) 

Jo 

The first order expansion of Eq. (10) is then given by 


Ps (t) =e tVl p s (t), ( 22 ) 

where 

V 1 =T(l + n)G- + TnG+, (23) 

G± = D[J±] — A± — B±. (24) 

For evolution under the first order average dissipator 
D i in Eq. (23) the expectation value of the ensemble 
magnetization is given by (See Appendix B) 

, N (l — e“*/ Tl ) 

(J z (t)) = e -t / Tl (J,(0))- 1 2 + 4 _ ’ ■ (25) 

where the relaxation time constant is given by 


T(1 + 2 n )' 


Q(t) = p(ir), G(t) = V cc (ir). 


(26) 
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This is an exponential relaxation process to an equilib¬ 
rium state with magnetization 

<•««• =-2T5- (27) 

Thus in the ideal cooling limit (n = 0) this is a T\ process 
to the ground state of the spin ensemble with magneti¬ 
zation ( J z ) eq = —TV/2, and Ti = T _1 . 


B. Simulation 

We numerically compared evolution under T>\ to both 
the cavity cooling with local deplrasing spin master equa¬ 
tion in Eq. (10) (Fig. 1), and to the full spin-cavity evo¬ 
lution under the Tavis-Cummings master equation with 
cavity dissipation and local spin dephasing (Fig. 2): 

j t p(t) = -i [. H TC ,p(t )] + {V c + V T2 )p(t). (28) 

All simulations were done using the QuantumUtils for 
Mathematica package [40] in the superoperator represen¬ 
tation for the system dynamics [41], using an TV-spin to¬ 
tally symmetric subspace matrix representation of the 
SU(4) algebra. In Fig. 1 we simulated cavity cooling with 
local dephasing for TV = 10 spins, and TV = 100 spins with 
dephasing rates of 7 = ATVT with A = 0,0.1,1,10. In 
both cases we see that for 7 = 0 we have cavity cooling 
alone and population is trapped. For 7 > 0 we break the 
SU(2) symmetry and achieve cooling to to the ground 
state, at a rate that increases with 7. For 7 = 10TVT 
we have good agreement with the first order average dis- 
sipator expansion V\. This is expected as Eq. (25) is 
only a valid approximation of the true dynamics when 
the dephasing strength 7 is sufficiently strong to disre¬ 
gard higher order terms in the Magnus expansion where 
terms that do not commute with the dissipator rapidly 
average to zero. In practice this corresponds to requiring 
7 > riV. The parameter C = rJV /7 is also called the co- 
operativity of the spin ensemble, and hence the condition 
for the validity of lowest order Magnus approximation is 
that C < 1. For example, in a typical X-band pulsed 
ESR setup reasonable values for cavity dissipation and 
single spin coupling are n/2ir = 1 MHz and g/2ir = 1 
Hz [18], and the Markovian condition for cavity cooling 
is satisfied for sample sizes TV < 10 12 [15]. For these 
parameters a single spin dephasing rate of 7 > 4 MHz 
would satisfy C < 1. We note that the relaxation curve 
is bi-exponential for C > 1. Once the initial collective 
cooling at rate T is saturated the long time rate is limited 
by the local dephasing rate. 

In Fig. 2 we simulated for TV = 10 spins and a cav¬ 
ity truncated to 4 levels with values of the spin-cavity 
coupling of g = 100, cavity dissipation rate k = 4 g 2 = 
4 x 10 4 , and spin dephasing rate 7 = A TV with A = 
0, 0.1,1,10,10 2 ,10 3 ,10 4 ,10 5 . The values of g and k were 
chosen to satisfy the Markovian condition n g\J TV for 


TV = 100, while giving an effective spin cavity cooling 
rate of T = 4g 2 /k = 1. In addition, the strong cavity 
dissipation rate allows us to truncate the cavity to low 
dimension. We find that the spin-cavity master equa¬ 
tion is in agreement with the spin cavity cooling master 
equation for dephasing rates up to lOTVr, however as the 
dephasing rate increases beyond the collective cavity dis¬ 
sipation rate, the cooling rate begins to slow down. Fol¬ 
lowing the derivation in [15], this can be incorporated by 
a Lorentzian cavity cooling rate T = 4 g 2 n/(K 2 + 4A 2 ), 
where A is a parameter that depends the the physical 
mechanism that gives rise to the local dephasing param¬ 
eter 7 . For example, in the simplest case of TV = 2, this 
is the expression where the dephasing mechanism arrises 
due to the spins being tuned to ±A/2 away from reso¬ 
nance with the cavity. 

The reduction in cooling rate in the strong dephas¬ 
ing regime is because the master equation in Eq. (10) is 
no longer an honest description of the dynamics: If we 
consider the full spin-cavity master equation in Eq. (28), 
then as the dephasing rate increases past the collective 
coupling strength this will suppress the TC exchange in¬ 
teraction for each individual spin with the cavity. To 
see this we note that in the dissipative interaction frame 
of T>t 0 , the spin cavity master equation in Eq. (28) has 
no time dependent piece (See Appendix C). Hence if we 
make an imaginary-time transformation and perform a 
first order Magnus expansion as done in Appendix B the 
exchange interaction will be completely averaged out to 
zero. Thus we find in terms of the full spin-cavity inter¬ 
action the local-XT cooling interaction enters at higher- 
order in the Magnus expansion. 


III. DISCUSSION 

The method we have described shows that cavity cool¬ 
ing techniques can be used to drive an ensemble system 
into the ground state by coupling to a cavity in the pres¬ 
ence of a local dephasing on the spin system. This tech¬ 
niques could prove useful for initialization an ensemble 
spin system in a highly pure state by short circuiting its 
thermal relaxation. There remains an opportunity for 
yet more efficient methods and the approach we have in¬ 
troduced for computing the cooling rates could be used 
more broadly. In addition the method for dissipative per¬ 
turbation theory that we developed to solve the cooling 
master equation should be useful for other systems where 
a dissipative term is dominant. 

Dephasing is always present in a real physical system 
and this proposal uses it as a resource for dissipative state 
engineering. Depending on the physical system there are 
many possible mechanisms that give rise to dephasing. 
For spins systems it may arises due to inhomogeneous 
static fields across the ensemble, which could be engi¬ 
neered by introducing gradient fields. For systems such 
as defects in silicon and nitrogen vacancy centres in di¬ 
amond it could arise due to hyperfine coupling between 
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N=10 N=100 



FIG. 1. Simulations of the magnetization expectation value (J z ( t )} for a maximally mixed initial state of an ensemble of N = 10 
spins (left), and N = 100 spins (right) for master equation described by first order average dissipator in Eq. (23) (dotted black 
line), and of the cavity cooling with local dephasing master equation in Eq. (10) for 7 = A1VT, and A = 0, 0.1,1,10. 



FIG. 2. Simulations of the magnetization expectation value ( J z (t )) for the full spin-cavity master equation in Eq. (5) with 
the addition of a local dephasing dissipator for N = 10 spins in a maximally mixed initial state, and a cavity truncated to 4 
levels initialized in the ground state. Evolution under the 1st order average dissipator is shown as the dotted black line in both 
figures. 


each electron and a local nuclear spin or spin bath, and 
the strength of these interactions can be engineered to a 
degree using decoupling sequences. 

Our calculation considered the case where the spin- 
ensemble was on-resonance with the cavity and the re¬ 
sulting cavity cooling dissipator can be thought of in 
terms of the Purcell effect. One could also consider the 
side-band cooling approach as detailed in [15] by intro¬ 
ducing a drive term on the spins to target a side-band of 
the resonator. In that case the magnetization of the spins 
under cooling accumulates in the J x basis (for a J x drive 
term) rather than the J z basis of the static field. In this 
situation the dephasing must also happen in the J x basis 
to achieve cooling to the J x ground state. In practice 
this could be engineered by using a gated protocol where 
a single cooling step consists of: side-band cavity cooling 
for a time i cc , applying a collective rotation swapping the 


J x and J z eigenstates, dephasing for time tr 2 in the J~ 
basis, then applying the inverse collective rotation to ro¬ 
tate back to the J x basis. This cooling step can then be 
repeated to form a discretized cavity cooling cycle with 
dephasing to reach the true ground state. 
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Appendix A: SU(4) Representation of Collective 
Spin Dissipators 
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We now outline how to represent the dissipators for an 
ensemble of N two-dimensional subsystems in terms of 
the generators of SU(4). Rather than explicitly derive 
this in terms of irreducible representations, we simply 
give a constructive method for representing each of the 
generators of SU(4) in terms of superoperators acting 
on the N subsystems. This approach has been used by 
[36, 37] and we use their notation for the SU(4) algebra 
here. 

There are 15 generators for SU(4), which each be¬ 
long to one of 6 SU(2) subalgebras. Let O = 
{Q, £, A4 ,MMi V} be the set of subalgebra operators. 
For O £ O we have 


[ 0 +, 0 -] = 20 3 , [ 0 3 , 0 ± }=± 0 ± . 


Also the pairs of operators (Q, £), (Ad, A/ - ), and {U,V) 
each commute: 



M+ 

Ad 

Ms 

A/+ 

M- 

A/3 

Q+ 

0 

-v+ 

-|2+ 

0 

U+ 

-12+ 

Q- 

V_ 

0 

|2- 

-U- 

0 

IS- 

Qs 

\M+ 


0 

W+ 

-W- 

0 

£+ 

0 

U- 

-|£+ 

0 

-V- 

-|£ + 

£_ 

-U+ 

0 

1 ^- 

V+ 

0 

|s- 

£3 

\M + 


0 

W+ 

-W- 

0 



U+ 

U- 

w 3 

V+ 

V- 

V 3 

Q+ 

0 

-A/+ 

-|s+ 

0 

A 4 + 

-|2+ 

Q- 

M- 

0 

IS- 

-Ad_ 

0 

|S- 

Qs 

\u+ 

-|w- 

0 

|V+ 

-|v- 

0 

£+ 

-M+ 

0 

§£+ 

A/+ 

0 

|£+ 

£_ 

0 

At- 

-§ E - 

0 

-A/L 

-|S- 

£3 

1 

toll-' 

£ 

+ 


0 

-|V+ 

|v- 

0 



U+ 


Us 

V+ 

V- 

v 3 

M+ 

0 

-£+ 

-\M+ 

Q+ 

0 

\M + 

M- 

£_ 

0 

1 M~ 

0 

-fi- 

~\M- 

Ms 


- 1 ^- 

0 

-|V+ 

|v- 

0 

A/+ 

-Q+ 

0 

1A/V 

0 

£+ 

— |A/+ 

A/L 

0 

Q- 

-W- 

-£_ 

0 

|A/"- 

M 3 

-|w+ 

W- 

0 

|v+ 

-|v- 

0 


[Qa^p] = [M a Mp\ = [U a ,Vfi\ =0 Va,/3e {±,3}. 


The remainder of the SU(4) commutation relations are 
shown in Table 1. Note that only 15 of these operators are 
linearly independent. In particular A/3 = Q 3 + £3 — .M 3 , 
U 3 = Ad 3 — £3 and V3 = Q 3 — Ad 3. 

The SU(4) generators may be given an explicit matrix 
representation in terms of superoperators acting on the 
vectorized density matrices of N two-dimensional sub¬ 
systems. This is analogous to how the spin operators, or 
SU(2) algebra, can be given a (2 J + l)-dimensional ma¬ 
trix representation in terms of the spin operators acting 
on a spin-J particle. Recall that for a quantum system 
with Hilbert space X = C d , density matrices are square 
matrices corresponding to linear maps p : X —t X. We 
can vectorize density matrices by stacking the columns 
together to form a column vector | p)) £ X®X. Superop¬ 
erators are then the linear operators which act on vector¬ 
ized density matrices: S\p)) = \ p')}. A linear map £(p) 
given by £ (p ) = ApB ! can be written as a superoperator 
as Ss = B*®A where * denotes complex conjugation [41]. 
A superoperator acting on N subsystems can be written 
as Sg = ^2jB^* g> A^\ and in particular the SU(4) 
operators may be represented in terms of the single spin 
Pauli operators a±\E±^ = |(1 + crl"^), and collec¬ 
tive spin operators J z = 02 ^/ 2 , J± = a ± ^ as 


TABLE I. Commutation relations for SU(4) algebra.The table 
entry is the value of the commutator [O tow , O co \] where O row 
and O co i are the corresponding operators in the same row and 
column of the table labels respectively. 


follows: 

N 

Q ± = £ , 

3 =1 
N 

s ± = E( 4 j) ® 4 J) )’ 

3 =1 
N 

M ± = J2 ( E + ® a ±) > 

3 =1 
N 

A4 = Y1 ( E - ] ® a ±) > 

3 =1 
N 

u± = E { a ± ® E +) > 

3 =1 
N 

v ± = E ( ct ± } ® E - ] ), 

j =1 


q 3 = J (1 ® J z + J z ® 1), 

£3 = * (1 <g> J z - J z <g> 1), 

1 N 

• m 3 = ^E(4 j) ®^ ) ) > 

3= 1 

1 n 

n, = 1 £ ( £ -’ ® d J) ). 

3=1 
* N 

u, = \ £ (d» ® e?) , 

3=1 

1 N 

v 3 = l £ (<*> ® e 0) ) ■ 

J=1 


The utility of this approach is that we can express 
many useful open system dynamics of collective spins in 
this representation, and it allows us to analytically com- 
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pute certain properties that may be difficult otherwise. 
In particular local T 2 and T\ process dissipators may be 
expressed as 

N 

Y jD [oW/2)=Mz- l -Qz-hlz-^l, (Al) 

3 =1 

N 

Y D W±] = 22± + 223 - NI, (A2) 

j=i 

where I = t is the SU(4) identity operator. Collective 
T 2 and T\ process may be expressed as 

D[J Z ] = - 2£§, (A3) 

D[J±] = (U± + V±)(M± + Af±) 

-\(Ut+V t )(U ± +V ± ) 

—-(Adip + A/) f )(AI± + A/±). (A4) 


Appendix B: Derivation of Local Ti Dissipation 

Consider a system with a cavity cooling dissipator, and 
a local T 2 dissipator. 

T> cc = r(l + n)D[J_\ + TnD[J + ], (Bl) 

N 

Vt 2 = 

3 =1 

The superoperators for these dissipators are given in 
terms of SU(4) generators by Eq. (A4) and (Al) respec¬ 
tively. We now consider the effective cavity cooling dis¬ 
sipator in the interaction frame of the T 2 dissipator: 


T> cc (t) = T{1 + n)D[J_]{t) + TnD[J + ](t), (B3) 

where 

D[J±](t) = e® T ^[4]e _ ® T 2. (B4) 

We may expand this using the BCH exansion: 

00 j-k 

D[J±}(t) = Y^ k [D[J ± }], (B5) 

k =0 

where Ck [7?[T±]] are nested commutator terms with 

C 0 [D[J±\] = D[J ± ], (B6a) 

Ci[D[J ± ]] = [Vt 2 ,D[J ± ]], (B6b) 

C k [D[J±]] = (B6c) 

We may compute the commutator terms of the BCH ex¬ 


pansion using the SU(4) algebra. To begin, we have 
[X>T 2 , M±\ = ±^Ad±, [Dt 2 , Af±] = (B7a) 

[Dt 2 , U±] =±\U ±1 [D Ti , V±] = t|v±. (B7b) 



Hence 


Vt 2 , D[J ± ] j = - \\Vt 2 , (G± + V±)(Ad± +A4)] 

+ [Dt 2 , (M±+Ar±){U±+V±)] 
- [V Ta , (G T +V T )(G± + V±)] 


= ±- 


- [Dt 2 , (M T +A7 T )(Ad ± +A4)]) 
(M±U± + U±M± - V T U ± - Af T M±) 


-(A4V± + V±A/± - G T V± - M t AT±) 


(B8) 


Define superoperators 
1 


A± = ~{M±U ± +U±M± - V T U± -AT t M±), (B9) 
B± = \ (A4 V± + V±A4 - U T V± - M T Af±), (BIO) 


then we may write 

Vt 2 , D[J±] = ±yA± =p , yB±. 


(Bll) 


If we then take the commutator of A± and B± with the 
T 2 dissipator we find 

[D T2 ,A±]= ± 7 A±, [T>t 2 ,B±\= =f 7 B±. (B12) 

Hence for k > 1 we have that the nested commutator 
terms are given by 

C k [D[J ± }] = (±7 ) k A ± + {Tl) k B ± . (B13) 

Thus the interaction frame dissipator terms are given by 

" t k 

k ! ( 


D[j±m = Y^c k [D[j ± }] 

k =0 

= D[J ± ] — A± — B± 


E 


(±7 ty 


A 4 


E 


(T7 tf 


B± 


k\ ^ ^ k\ 

k=0 k—0 

= D[J±] + (e ±7t - 1)A± + (e ** - 1 )B±. 

(B14) 

If we perform a change of variables to imaginary time 
the interaction frame operator G(t) = D[J±\(ir) is pe¬ 
riodic with period T = 2'rt/ r )~ 1 . The first order Magnus 
term over this period is the time-independent piece of the 
dissipator: 


V 1 =G 1 = T(1 + n) G- + TnG + , 


(B15) 


where 


G± = D[J±] — A± — B± 

= l (u±M± +Af±U± + Ad±V± + V±Ad±) 

— — :W± + VipV± + AdzpAdi + AfzpAf±'j. 

(B16) 














1. J z Expectation Value 


and 


We are interested in computing the evolution of the 
expectation value for the J z operator for an arbitrary 
initial state for dynamics described by the first order av¬ 
erage dissipator in Eq. (B15). We do this by solving 

(j z (t)) = « JMt )» = ((JzW))) 


N 

[U+M+ + M+U+) | J z )) = Y | E + o-+ Jz ° { - )) 
i,j =1 


N 


i,j =1 

= | (f 1 + J z ) (J, - 1))) 


where we have made use of the relations 


where J z {t) = J z (t) as [J z ,T>t 2 ] = 0. The adjoint dissi¬ 
pator is given by 


v\ = G[ = T(1 + n)Gl + rnG f + , (B17) 


where 


g \=1 


(u t a r T + u T u T + v t + 

— — + VzpV± + AdzpA 4± + AfzpAf±^j , 

(B18) 


and we have used Ad± = Ad=p and similarly for 

Jf±M± ,V±. To compute the terms of G^| J~)) we have 
that 


AT 


+ A/’zpWrp)! J z )) — 

bi=t 


N 

= E 

*d=i 


4 j) 4 i)j ^- )f7 ± } 


4 ) 4 J 4°± ) £- ) 




4 )</ * £7 + )<r ± )£7 - ) 


Hence 


N 

{U-U- +AT-U-)\ J z )) = Y | ^ 0+ 5 s - ) }) 

*»j=l 
N 

= E 

*b=i 


J* 


4 } ^ 


JzO-^4^ 


N 

E| (Jz + i)^)) 


j=i 


= l(j,+i)(fi- j z ))>, 


E±E± = E±, 

E±Ez p = 0, 

E±a± = cr±Ezf = cr± 

Ez p<r± = cr±E± = 0. 

Thus we have 

J z » = | (f 1 =F Jz) (Jz ± 1 ))). 

Similarly one can show 

(V T Af T +Af T V T )|Jz)) = {(flT Jz)(Jz±l)»- 

For the other terms we have 
N 

U T U ± \J Z )) = Y | E^E^J^a^)) 

i,j =1 

N 

bt =1 

and hence 

N 

UM +1 J 2 )) = E | S + } ^+^-^+ ) )) 

= 0 , 

N 

U+U-\J Z ))= Y | Jz^^o- 5 )) 

*»j = l 
N 

= £|/.£?>)) 
i=i 

= |jz(fl + Jz))). 

Similarly 

Ad_AJ+| J 2 )) = 0 , 

Ad+Ad_|Jz)) = |Jz(fl + Jz))), 
V_V+|Jz)) = |Jz(fl- Jz)», 

V+V_|Jz)) = 0, 

V-V+l Jz)) = | Jz(fl- Jz))), 

V+V-l Jz)) = 0. 
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Thus we have 


G\\J Z )) = | (f 1 + J z ) Jz)) ± | (f 1 T Jz))) 


- | (f 1 + Jz) Jz 
= -|J,»±^|l}>. 


(B19) 


T^t 


Next we need to evaluate G ± | 1)). We have that 


A*+| 1 » = 0 , 

U+ | 1 »= 0 , 

W + |1)) = |J + )), 

V + |1)) = |J_}), 


^_|1» = |J_», 
W_|1)) = |J+)), 
AA_|1» = 0, 
V-|1» = 0, 


Hence the effective dynamics are described by an expo¬ 
nential decay process with decay rate 


Ti = 


1 

T(1 + 2 n) 


(B24) 


to an equilibrium state with magnetization 


{Jz) 


eq — 


N 

2 + 4n’ 


(B25) 


where in the ideal cooling limit of n = 0, this is a Tf 
process to the ground state of the spin ensemble. 


and so 

G^|l)) = i(AL|J+)) + V_|J_))) 

-i(A/-_|J + » + V-|J_») 

= 0 , 

Gl\t))= l -(M + \J-))+U + \J + ))) 

- 1 -[m + \j.))+u + \j + ))) 

= 0 . 


Finally we may put this all together to obtain 


7=4, 




G | J z » = F (Jl + n)Gl +nG \J | J z )) 

= r(l + n)(-|J,))-y|l))) 

N, 


+ Tn( —| J z )) + —\ 1») 


= -r(l + 2 n) 

and in general for k > 1 


Jz)) 


N 


2 + 4 n 


| 1 » 


k 


G \J Z )) = (-1) r (1 + 2n) 


Jz)) + 


(B20) 




(B21) 


Appendix C: Strong Dephasing Limit 

In this appendix we outline a model for the cooling 
results in the presence of strong dephasing in Sec. IIB. 
The cooling dynamics in Fig. 2 shows the spin-ensemble 
magnetization for evolution under the the full spin-cavity 
master equation with local dephasing: 

i\p { t))) = (s TG +v c +v T 2 )\ P m, (ci) 

where p{t) is the density matrix for the joint spin-cavity 
system and <Stc is the superoperator for the TC interac¬ 
tion 

I p(t)')} = 5tc| p(t))) -<=> p(t)' = -i[H TC ,p(t)\. (C2) 

As described in Sec IIB, as the dephasing rate in¬ 
creases beyond the collective cavity dissipation rate the 
cooling rate begins to slow down. This is because in the 
strong dephasing regime the spin-cavity interaction term 
6>tc is suppressed. If we consider the Magnus expansion 
of <Stc in the interaction frame defined by T>t 2 we find 
that there is no secular piece that commutes with the de- 
plrasing interaction frame. Let S± ' be the superoperator 

defined by «S±^|p s )) ^ —i[cr±\ps]. In column stacking 
convention this is given by 


Hence we have 

oo j-k 

k—0 


= e -tr(l+2n)|j^ 


'n( 1-e 


—t r(l+2n 




2 + 4n 


| 1 »- 

(B22) 


Thus for an arbitrary initial state p the expectation value 
of J z under this evolution is given by 


(J z {t)} = e~ tr(1+2n \j z ( 0)) 

_ (l _ e -tr(l+2n) 


N 


2 + 4 n 


(B23) 


S £ j) = -i (a® ® l (i) - . (C3) 

If we move into the interaction frame of r )'D[a z ^] we have 

S(t) = -i (e 2 * 7 ® 1^) - e" 2t7 1^) (8> a ( J^ (C4) 

and so all terms of <5>tc(£) have time-dependence of e ±2t7 . 
We can see here that if we make an imaginary-time trans¬ 
formation and perform a first order Magnus expansion as 
done in Appendix B that the exchange interaction will 
be completely averaged out to zero. Hence in terms of 
the full spin-cavity interaction the local-T) cooling inter¬ 
action enters at higher-order in the Magnus expansion. 
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